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Possible candidates for the Mg2B30i: nanostructures observed in bulk of poly crystalline MgB2 
Pi have been studied using a combination of Z-contrast imaging, electron energy loss spectroscopy 
(EELS) and first-principles calculations. The electronic structures, phonon modes, and electron 
phonon coupling parameters are calculated for two oxygen-ordered MgB2 compounds of composition 
Mg2B30 and Mg2B302, and compared with those of MgB2. We find that the density of states for 
both Mg2B30i structures show very good agreement with EELS, indicating that they are excellent 
candidates to explain the observed coherent oxygen precipitates. Incorporation of oxygen reduces 
the transition temperature and gives calculated Tc values of 18.3 K and 1.6 K for Mg2B30 and 
Mg2B302, respectively. 



The discovery of superconductivity in MgB2 with a 
transition temperature of 40 K has attracted the at- 
tention of the scientific community for two main reasons: 
the technological applications of this material and the 
new insights that such a simple structure could bring to 
superconductivity theory. Efforts to improve the super- 
conductivity properties of MgB2 have included doping 
with elements such as Y 0, Zr 0, C [|, Al [|, Cu 
and Ag 0, but the results varied for different groups. 
A simple explanation for this behavior is that the trans- 
port properties of MgB2 have a strong dependence on 
the sample preparation conditions, which makes it hard 
to obtain unambiguous results even for undoped samples 
0. In fact, there have been reports from different groups 
showing that MgB2 samples are not a single phase mate- 
rial, but rather a rich collection of different phases includ- 
ing MgB2, MgO, ByO,, Mg,B^, Mg.B^O, 
Oxygen, unlike some other elements, is present in MgB2 
as an unintentional impurity due to its high reactivity 
with MgB2. Typically, oxygen rich precipitates form in 
the bulk of MgB2 where they contribute to the overall 
transport properties through flux pinning [llj. Moreover, 
oxygen can segregate at grain boundaries of MgB2 where, 
unlike the high— Tc materials, it contributes significantly 
to flux pinning, increasing the overall Jc as the grain size 



decreases MgB?. has already been used in making 

wires 0, tapes 0, and Josephson junctions J^j, and 
in order to improve the properties of such devices, the 
effect that oxygen has on the properties of MgB2 needs 
to be well understood. Surprisingly, the importance of 
oxygen in the development of MgB2 devices has been 
overlooked, and detailed first-principles studies of oxygen 
ordering and segregation have not yet been performed. 

It was previously reported that oxygen impurities can 
increase the upper critical fields and critical current but 
lower Tc [l^ . Liao et. al. found that through a slow 
cooling process, the oxygen in bulk MgB2 can segregate 
to form nanometer-sized coherent Mg(B,0) in the MgB2 
matrix actin g as effective flux pinning centers without de- 
creasing Tc . Klie et. al. previously reported oxygen 
precipitates in bulk MgB2, where oxygen was incorpo- 
rated every second boron plane of the MgB2 structure, 
found by atomic resolution electron energy loss spec- 
troscopy (EELS) studies using a scanning transmission 
electron microscope Ij. Figure 1 shows an atomic res- 
olution Z-contrast image of one of these oxygen precipi- 
tates taken from bulk MgB2 along the [100] orientation. 
The increase of contrast every second column in Fig. 1 
is due to the presence of oxygen 0|. The size of these 
precipitates was reported to be about 20-100 nm 0. 
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FIG. 1: Z-contrast image of a coherent oxygen precipitate in 
the bulk of MgB2 [010]. The bright spots are Mg atoms. The 
increase of contrast every second column is due to presence 
of oxygen every second boron plane. White arrows highlight 
the atomic columns where oxygen has precipitated. 



Motivated by these experimental observations, in this 
work we discuss the effects of oxygen-ordering, with dif- 
ferent stoichiometrics, on the electronic structure and 
superconductivity properties of MgB2 via first-principles 
calculations. As possible candidates for the atomic struc- 
tures of coherent oxygen precipitates in bulk MgB2, we 
propose two different structures of composition, Mg2B30 
and Mg2B3 02. We find that incorporation of oxygen in 
the form of Mg2B302; structures decreases the DOS at 
the Fermi energy and the electron phonon coupling com- 
pared with MgB2. From total energy calculations we 
find that Mg2B30 is less stable compared to Mg2B302 
but has a higher Tc- It is important to notice that our 
calculations show that oxygen impurities are likely to be 
present everywhere in MgB2 structures. Therefore, con- 
trolling the density and type of these impurities will fi- 
nally be the limiting factor in improving the transport 
properties of MgB2. 

Our atomic and electronic structure calculations were 
performed using the ab initio pseudopotential plane wave 
method [l^ within the generalized gradient approxima- 
tion (GGA) 13, and local density approximation (LDA) 
for comparison. In order to compare the changes due to 
oxygen in MgB2, we doubled the unit cell of MgB2 in 
the c-direction. Brillouin-zone integrations were carried 
out with 15x15x6 Monkhorst-Pack grids for Mg2B4 and 
Mg2B302;. An energy cutoff of 400 eV was used for all 
the calculations. The lattice volumes, shapes, and all 
atoms were allowed to relax for the new structures using 
a force criterion of 10^^ eV/A. The zone-center phonons 
and electron phonon couplings were calculated using den- 
sity functional perturbation theory j2oj . 



FIG. 2: (color). From left to right: schematics of Mg2B4, 
Mg2B30 and Mg2B302 crystal structures. The lattice in the 
c-direction stretches for Mg2B302 structures with respect to 
Mg2B4, where the addition of oxygen increases the distance 
between the Mg planes where oxygen was added. Mg, B, and 
O atoms are shown in gray, yellow, and red respectively. 

The relaxed geometries for structures with partial sub- 
stitution of boron atoms with one and two oxygen atoms, 
Mg2B30 and Mg2B302, as well as Mg2B4, are shown 
in Fig. 2. Other oxygen concentrations were tried as 
well (i.e. Mg2B303, where one additional oxygen was 
inserted in the boron plane), but these resulted in a sig- 
nificant distortion of the MgB2-like symmetry and did 
not reproduce features observed in the coherent oxygen 
precipitates. We found that the incorporation of oxygen 
changes the well-defined B-B planes by the formation of 
BOa; units in the Mg2B302; structures. Both Mg2B30a: 
structures belong to the space group P3ml (Cg^,, No. 
156). The relevant structural parameters for the three 
structures considered are shown in Table Introduc- 
tion of oxygen considerably increases the distance be- 
tween neighboring Mg planes, but decreases the distance 
between the other Mg planes, which are above and be- 
low the boron layer left intact. As a result, the lattice 
stretches in the c-direction compared to Mg2B4. Such be- 
havior was reported by Eom et. al. where thin films rich 
in oxygen showed a higher c lattice parameter suggest- 
ing that boron can be substituted by oxygen 116|. The 
calculated ratios of the distances between the Mg-BO^,- 
Mg and Mg-BB-Mg layers are 1.22 and 1.49 for the the- 
oretical Mg2B30 and Mg2B302 structures, respectively, 
while the estimate from the experimental image is 1.14 ± 
0.05. The overestimated theoretical values for the ratios 
compared with the experimental can be explained in two 
different ways: (i) pressure effects from the surrounding 
bulk into the precipitate, which shrinks the lattice pa- 
rameter in the c-direction, and (ii) the concentration of 
oxygen in the BOa; planes is slightly lower than the ones 
calculated here, where x < 1. 

We also investigated the phase stability of the 
Mg2B30a; compounds relative to atomic constituents and 
bulk MgB2. Both Mg2B302; structures are stable against 
phase separation into their atomic constitutes. How- 
ever, the magnitudes of the binding energies per atom 
for Mg2B30 are found to be 0.16 eV and 0.63 eV smaller 
compared to that of MgB2 calculated using atomic oxy- 
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TABLE I: Summary of properties for MgB2, Mg2B4, Mg2B30 and Mg2B302 crystal structures. The Mg2B30i: structures 
belong to the space group P3ml (Cj^,, No. 156). Lattice parameters are in A, N{Ep) in states/eV, frequencies (lo) in meV. 
The parenthesis next to the frequencies indicates degeneracy of those modes. Tc is given in K. The lattice parameters are 
shown for GGA (LDA) approximations. The phonon modes, electron phonon coupling, and the critical temperatures were all 
calculated using the values from the GGA calculations. 



MgB2 


Mg2B4 




Mg 


2B3O 


Mg 


2B3O2 


a^s.ose" 


a=3.046(3.086) 




a=3.067(3.016) 


a=3.143(3.016) 


c=3.524'" 


c=7.054(7.048) 




c=7.701(7.545) 


0=8.633(8.231) 




Ai'(^SF;=1.35(1.38) 


N(Ef)= 


=0.85(0.87) 


N(Ef)= 


=1.05(1.99) 


LUp Ap 




Ap 


LOp 


\ 

/\p 


LUp 


\ 

y\p 


87.1(l)Bi9 


26.3(2) 





22.5(2) 


0.027 


10.0(2) 


0.047 


74.5(2)E29 0.907 


27.0(2) 


0.06 


28.4(2) 


0.029 


16.8(2) 


0.074 


49.8(1)A2„ 


39.0(1) 


0.08 


39.2(1) 


0.109 


21.0(2) 





40.7(l)Ei„ 


39.7(2)Ei„ 





39.3(2) 


0.010 


25.0(2) 


0.003 




47.6(1) 





48.9(2) 


0.026 


33.6(1) 







49.3(1)A2„ 





51.0(1) 


0.019 


36.5(1) 


0.124 




68.0(2) 


0.51 


53.0(1) 


0.005 


41.2(2) 


0.002 




73.4(2)E2g 


0.46 


72.1(2) 


0.506 


55.4(1) 







79.2(1) 


0.03 


90.4(1) 


0.002 


83.1(1) 







86.3(l)Bi9 





161.8(1) 


0.088 5 


89.121) 


0.153 












105.3(1) 


0.131 












215.8(1) 





A = 0.907 


A = 1.14 




A = 


0.821 


A = 


= 0.534 


(Jlog = 74.5 


U}log = 64.5 






= 65.6 


LOlog 


= 48.9 


Tc = 39.4 


Tc = 40.73 




Tc -- 


= 18.31 


Tc 


= 1.62 



"Taken from Ref. 2] 

''Values for N(Ep ), u), A and Tc (with ^j* = 0.15) are taken from Ref. [g^ 



gen, and molecular oxygen, respectively. On the other 
hand, the magnitudes of the binding energies per atom 
for Mg2B302 are 0.11 eV higher, and 0.71 eV lower com- 
pared to that of MgB2 calculated using atomic oxygen, 
and molecular oxygen, respectively. These calculations 
suggest that (i) the theoretical Mg2B30a; structures are 
not stable against phase separation into MgB2 and molec- 
ular oxygen, and (ii) Mg2B302 is stable with respect to 
phase separation into MgB2 and atomic oxygen. This 
does not, however, mean that these structures have to be 
ruled out as possible candidates for coherent oxygen pre- 
cipitates, since they can still exist as metastable phases 
in bulk MgB2. 

Figure 3 shows the calculated unoccupied boron DOS 
for the theoretical Mg2B30a; nanostructures compared 
with the experimental boron K-edge spectra taken from 
the B-B plane (darker columns in Fig. 1) and BO^; plane 
(brighter columns in Fig. 1) respectively. Each experi- 
mental spectrum is the sum of seven individual spectra 
with an acquisition time of 1 sec. and an energy reso- 
lution of 1.2 eV. The spectra have been summed up to 
increase the signal to noise ratio, and background sub- 
tracted and corrected for multiple-scattering contribu- 
tions To obtain a better comparison between theory 
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FIG. 3: (color). Experimental Boron K-edge spectra of co- 
herent oxygen precipitate (black) compared with calculated 
density of states of Mg2B30 (blue) and Mg2B302 (red) (a) 
B-B plane, and (b) BOx plane. Solid (dashed) lines have an 
energy resolution of 1.2 (0.4) eV. Experimental monochro- 
mated spectra (green) with an energy resolution of 0.3 eV 
are shown for comparison for (a) bulk MgB2, and (b) a pre- 
cipitate region with oxygen. Notice that the pre-peak in the 
Boron K-edge is completely resolved for (a) and (b) . Changes 
due to oxygen in the monocromated spectra (b) are clearly 
notorious, change of intensity in the pre-peak and the peak 
at 195 eV. 
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and experiment, we have calculated the boron p^y and 
Pz states. We then took into account the scattering mo- 
mentum transfer dependence of the electron microscope 
optical conditions and crystal orientation to weigh the 
boron p^y and Pz states according to the experimental 
setup 21]. Analysis of the spectra for energies in the 
range of 15-20 eV above the Fermi level was not per- 
formed, since multiple scattering and lifetime broadening 
limit the resolution of the spectra, which makes compar- 
ison between experiment and theory uncertain in that 
regime [2^. The theoretical spectra was obtained in the 
Z approximation, which has worked well for bulk MgB2 
[2]|. have been convoluted for an energy resolution of 1.2 
eV and 0.4 eV (shown for comparison purposes only). 
Figure 3(a) shows the experimental and theoretical spec- 
tra for the B-B plane, where both structures, Mg2B30 
and Mg2B302 (energy resolution of 1.2 eV), have good 
agreement with experiment. The pre-peak of the boron 
K-edge at 189 eV in the experiment is well reproduced 
in both theoretical structures. A peak at 192 eV, which 
is present in both theoretical spectra, is not resolved in 
the experimental spectra. This is possibly due to the low 
signal to noise ratio of the experimental spectra. Fig- 
ure 3(b) shows the spectra for the BOa, plane, where the 
enhancement of the peak at 195 eV, which was previ- 
ously attributed to oxygen P, 0, |23| , is very well pre- 
dicted by both theoretical structures. This peak has a 
lower intensity for the B-B plane spectra (Fig. 3(a)). Fi- 
nally, the pre-peak in the BO^; plane experimental spec- 
tra shows an enhancement in intensity compared with 
the B-B plane experimental spectra. Here, both theo- 
retical spectra present a remarkable agreement with the 
experiment, especially the Mg2B30 structure where the 
difference in intensity of the pre-peak between the B-B 
plane and the BOj, plane spectra is more pronounced. 
Additionally, monochromated spectra (energy resolution 
of 0.3 eV) from bulk MgB2 and a region containing oxy- 
gen are shown in Fig. 3(a) and Fig. 3(b), respectively. 
Notice that the pre-peak is completely resolved for Fig. 
3(a) and Fig. 3(b). Furthermore, the monochromated 
spectrum (Fig. 3(b)) presents an enhancement due to 
oxygen of the peak at 195 eV. 

To investigate the effect of oxygen incorporation by 
partial boron substitution on the superconducting prop- 
erties of MgB2, we first calculated the phonon modes 
and electron phonon coupling parameters for Mg2B4, 
Mg2B30, and Mg2B302. The results are shown in Ta- 
ble In order to obtain the transition temperatures Tc 
of the theoretical Mg2B30a; nanostructures we used the 
McMillian expression [2^ 



-1.04(1 + A) 



Tn = CXp , 

^ 1.2 ^ \x- ^l*il + o.62^l*x) 



(1) 



with Xp being the electron phonon coupling for each mode 
at r, and fi* represents the Coulomb repulsion interac- 
tion. We used the value of /i* = 0.19. In the case of 
MgB2, the Eliashbcrg function has a sharp peak around 
the E2g mode ,2^; therefore, it is sufficient to take the 
mode energy as ujiog to obtain the transition tem- 
perature Tc [23, 113 • In fact, using only this mode, not 
only Tc but also the isotope effect and pressure depen- 
dence of Tc can be successfully explained However, 
for the Mg2B30a; nanostructures, the E2g mode is no 
longer well defined and uJiog has to be averaged through- 
out all the phonon modes. We approximated uJiog as 

exp ^j^^n {ujp) Xp^ , and only calculated the 

phonons at the zone center F. We obtained transition 
temperatures Tc of 40.73 K, 18.31 K, and 1.62 K for 
Mg2B4, Mg2B30, and Mg2B302, respectively. 

In conclusion, we have proposed two different struc- 
tures, Mg2B30 and Mg2B302 as possible candidates for 
the Mg2B30a; nanostructures observed experimentally in 
bulk MgB2 1] . Both structures present an increase of the 
c lattice parameter with respect to Mg2B4, where the 
Mg2B30 structure is closer to the experimental value. 
The experimental spectra for the B-B plane and BOj, 
plane are well explained by the Mg2B30 and Mg2B302 
structures. Both structures are good candidates for co- 
herent oxygen precipitates, although their superconduc- 
tivity properties differ. Mg2B302 precipitates would be- 
have as a pinning center for temperatures higher than 
2 K, while Mg2B30 precipitates would start to behave 
as a pinning center for temperatures higher than 18 K. 
Our results show that oxygen impurities are likely to be 
present everywhere in MgB2 , highlighting the importance 
of oxygen as an agent in the creation of flux pinning cen- 
ters and as a factor of improvement upon the critical 
current density Jc- 
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where ujiog is the logarithmic average frequency, A is the 
total electron phonon coupling defined as A = ^ Xp, 

per 
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